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Relaxation of a non-equilibrium state in a disordered metal with a spin-dependent electron energy 
distribution is considered. The collision integral due to the electron-electron interaction is computed 
within the approximation of a two-quasiparticle scattering. We show that the spin-flip scattering 
processes with a small energy transfer may lead to the divergence of the collision integral for a 
quasi one-dimensional wire. This divergence is present only for a spin-dependent electron energy 
distribution which corresponds to the total electron spin magnetization M = and only for non-zero 
interaction in the triplet channel. In this case a non-perturbative treatment of the electron-electron 
interaction is needed to provide an effective infrared cut-off. 
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1. Introduction. The ground-breaking experiments 
by Pothier et al. of Ref. [l| have demonstrated that one 
can have a direct access to the non-equilibrium electron 
energy distribution function f{E) and through it to the 
inelastic collision integral /Ccoii(^) which enters the ki- 
netic equation^: 

dtf{E;x,t) 



D\/^f{E;x,t) = -JC,oii{E;x,t). (1) 



In turn, studying the collision integral gives one an 
important information on interaction and dynamics of 
quasiparticles in a dirty metal. In this way the predic- 
tions of the theory of electron interaction in disordered 
metals^'^ were checked^^^ and an unexpected strong sensi- 
tivity of the energy relaxation to the presence of magnetic 
impurities^ was established. 

The main idea of Ref. 1 was to use the sharp fea- 
tures in the energy dependence of the density of states 
(DoS) ^probe(-^) of ^ supcrconducting probe electrode, 
which enabled to extract f{E) by measuring the differ- 
ential conductance of the tunnel junction between the 
normal metal sample and the probe electrode. Recently, 
the same idea has been suggested? to create a non- 
equilibrium spin- dependent electron energy distribution 
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FIG. 1: The setup: a quasi-one dimensional disordered 
normal-metal wire connected through the weak tunnel junc- 
tions to the two superconducting leads. The non-equilibrium 
is created by applying the finite bias voltage U between the 
leads. 



fa{E) and thereby to obtain a spin-polarized current 
through the probe. The sketch of the experimental setup 
is shown in Fig. [H The sample in a form of a quasi-one 
dimensional disordered normal-metal wire is connected 
through the weak tunnel junctions to the two supercon- 
ducting leads. The non-equilibrium is created by apply- 
ing the finite bias voltage U between the leads. The spin 
dependence of fa{E) {a = ±1 for the spin projections 
I (i)) is caused by magnetic field applied to the supercon- 
ducting leads with the DoS v^^\e) = us{E^U/2±aEz) 
and iy^J'\E) = vs{E - U/2 + crEz) where ys{E) = 
[E/^/E^ - /S?\ and the Zeeman shift Ez = IJLbH 
is taken with the sign zb depending on whether the di- 
rections of the magnetic field in the right and left leads 
are parallel or anti-parallel. In the absence of relaxation 
f^{E) is given by^^^^: 



{E)fF{E - U/2) + v^J'\E)fF{E + U/2) 



{E) 



(E) 



(2) 

where friE) is the Fermi distribution function. The 
measured quantity is the differential conductance with 
respect to the probe bias Vprohe across the probe tunnel 
contact. The probe contact can act as a spin- analyzer 
provided an additional magnetic field is also applied to a 
superconducting probe electrode. 

There are two distinct cases schematically shown in 
Fig. [2^ and Fig. [2)3: (i) with parallel and (ii) with anti- 
parallel magnetic fields in the superconducting leads. In 
the former case a non-equilibrium state with a nonzero 
total spin polarization 



M = I dE[f^{E)-f^(E)\ 



(3) 



is created, while in the latter case = 0. The typical 
form of the difference /dif = f}{E) — fi{E) that follows 
from Eq. ([2]) is shown in Fig. [2] in both cases. 
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FIG. 3: Three possible processes, allowed by conservation 
laws. In the initial state the quasiparticles have: (i) the same 
spin projection which remain unchanged during the collision, 
(ii) opposite spin projections which do not change during the 
collision, (iii) opposite spin projections and the collision re- 
sults in a spin- flip of both quasiparticles. 



FIG. 2: Two distinct cases: (a) with parallel and (b) with 
anti-parallel magnetic fields in the superconducting leads. In 
(a) case a non-equilibrium state with a nonzero total spin 
polarization is created, while in (b) case = 0. The typical 
form of the difference /idif = (^) ~ (^) shown in both 
cases. 



The aim of this work is to consider the relaxation 
of such a spin-dependent distribution caused by the 
electron-electron interaction. For this we derive the col- 
lision integral /Ccoii in the approximation of the two- 
quasiparticle collisions in the case where both the spin- 
singlet and the spin-triplet channel of the interaction are 
present. The detailed derivation of the collision integral 
due to the electron-electron interaction has been recently 
carried out in Refs. However, the analysis has been 
limited to the case of the spin-independent distribution 
functions /(^), while we are going to focus on the relax- 
ation of the difference /dif The results of calculation 
of the collision integral for the spin-dependent distribu- 
tion function have been also recently reported in Ref. 10. 
However, the authors considered a limited class of distri- 
butions with a very particular spin-dependence equiva- 
lent to a shift in the energy (E) = fi{E-\- SE) . Such 
type of dependence does not hold e.g. in the experimen- 
tal setup of Figs, m [21 

The main qualitative result of our analysis is that there 
are three different contributions in the collision integral. 
Two of them are also present if only the singlet channel of 
the interaction is considered, with only their amplitudes 
depending on the triplet channel interaction constant F. 
The third contribution corresponding to the spin- flip pro- 
cess is only present when the triplet channel of the inter- 
action is taken into account. Its magnitude dependsi^ 
essentially on the conserving total spin Ai. However 
more importantly, it is singular for the non-equilibrium 
spin-dependent distribution with JH = which naturally 
arises in the experimental situation (ii) of Fig. [2)3. The 
existence of such a singularity which never occurs if only 
the singlet channel of the interaction is present, is the 
main qualitative result of this work. 

2. Three contributions to the coUision integral. 
For a generic two-quasiparticle collision in a disordered 
metal in the absence of spin-orbit interaction and mag- 



netic impurities two quantities are conserved: the to- 
tal energy S and the total spin J\4. The latter con- 
servation law allows only three possible processes (see 
Fig. [3]): (i) in the initial state the quasiparticles have the 
same spin projections (l/2)cr which remain unchanged 
during the collision, (ii) in the initial state the quasi- 
particles have opposite spin projections which do not 
change during the collision, (iii) in the initial state the 
quasiparticles have opposite spin projections and the 
collision results in a spin-flip of both quasiparticles. 
Each process corresponds to a certain term in the col- 
lision integral that contains combinations of the type 
r{E + CO) r\E') [1 - f-{E)] [1 - f-\E' + u)] - [1 - 
/"(^ + ij)] [1 - f "'(^0] /^''(^) ^'''{E' + uj) which for 
the processes (i)-(iii) take, respectively, the forms: 

= j dE' [-{I - K^eK,E+oj){K^E'+oj - K^E') 

+ {ha,E+uj-K,E){l-K,E'K,E'+u)], (4) 

4'^ = JdE' [-{l-K,EK,E+oo){h-a,E'+oo-h-,^E') 

+ (hcr^E+uj — ^cr,£;)(l - h-a^E'h-cr.E'+uj)] , (5) 

= j dE' [-{I - K,Eh-a,E+u;){h-a,E'+uj - K,E') 

+ {h-a,E^u - K,e){^ - K,E'h-a,E'^u)] , (6) 

where h^^E = 1 - 2/a(£^). 

The collision integral JCco\\{o'^ E) can be represented as 
follows: 

3 

/Ceou(cT,i?) =PJ^ K;{lo)I(/HE,u;), (7) 

where u is the DOS (per spin direction) of the normal- 
metal sample at the Fermi level. The quantities Kp{u;) 
describe the strength of relaxation due to the correspond- 
ing processes (the quantities corresponding to the singlet 
channel 2{^) — ^1,2 (^) does not depend on Ad and 
hence on cr). We obtained the following expressions for 
them valid in the limit pf^ ^ 1 {pf is the Fermi mo- 
mentum, I is the elastic scattering length) and for the 
diffusive quasiparticle dynamics: 

K - ^ ^^ + ^)'+ (2^^ (o^ 



3 



-y- 

q 



(2Dq2) 



2T2" 



(1 + F)2 



(9) 



(Dq2)2 



^3H = 



q 



(a;-crFA4)2 

2)2 



(10) 

In Eqs. ©-(do]) by F (F = -1 corresponds to the Stoner 
instability) we denoted the Fermi-hquid constant corre- 
sponding to the triplet channel of the electron-electron 
interaction. The summation over q can be replaced by in- 
tegration ^^q^ J in the limit co > Erh = D/L'^ 
{D is the diffusion coefficient, L is the length of the dis- 
ordered sample) which will be considered below. 

Note that at small F <C 1 we obtain up to the linear 
in F terms: Ks = and 



-E 



1± 



2F 



1 



(11) 



In this limit the spin-spin interaction in the triplet chan- 
nel results in only a small (and opposite in sign for the 
parallel and anti-parallel spins of the two quasiparticles 
in the initial state) change in the amplitudes of the pro- 
cesses (i) and (ii) which are dominated by the singlet 
channel of the electron-electron interaction. Note that 
under the restrictions on the form of the spin-dependence 
of hfj^E adopted in Ref. 10 the combinations /^^^ and 
la appeared to be identical. This is why the result of 
Ref. [2 contained only the combination {Ki + K2)I^^ 
and K^I^a^. 

3. Relaxation of a non-equilibrium distribution 
and the conservation laws. As has been already men- 
tioned, the form of the collision integral should be com- 
patible with the two conservation laws. The conservation 
of the total energy requires: 



dE E [/Ccon(T, E) + /Ccon(i, E)] = 0. (12) 
The conservation of the total spin polarization leads to: 



d£[/Ccon(T,-E)-/Ccona,-E)]=0. 



(13) 



One can check using Eqs. (|4|)- (p!Q|) that fulfillment of 
those two conservation laws is guaranteed by the struc- 
ture of I^\e^ u) and the following properties of the ker- 
nels Kp{u;): 

Ki^2{co) = Ki^2{-co), Kiiu) = K^%-uj). (14) 

What we would like to note here is that the full relaxation 
to equilibrium due to the electron-electron interaction is 
only possible if = 0. Indeed, the fixed solutions to the 
kinetic equation Eq. ([T|) which correspond to all combina- 
tions Eq. (|4])-(|6]) vanishing identically, are the Fermi dis- 
tribution functions /j^^(F) = /^(F ^ M/2). Any non- 
equilibrium distribution tends to relax to these fixed solu- 
tions. However only at = we have /i^^ (F) = /[^^ (F) 



which corresponds to the complete equilibrium. So we 
encounter for the first time with the special role of the 
= condition. 
4. Collision integral for a quasi- Id wire and the 
infrared catastrophe at = 0. For the quasi- Id 
experimental geometry of Fig. [H Eqs. ([8|)- (p!Q|) can be 
straightforwardly evaluated: 



i^l,2(c^) = 



1 



^3(^) = 



85V'2F(1 + F) 



13/2' 



(15) 



(16) 



s {^\uo-(jFM\ + viM vTTf) 
1 

[|c^ - (jFM\ + \uj\{l + F)] V2F(1 + F) ' 
where S is the cross-section area of the quasi- Id wire. 



1 



4F(1 + F) 



(1 + Vr+F)(2 + F) 



(17) 



C2 = 1, Cs = F^/2 and the legitimate values of F are 
F > -1. 

A remarkable property^^ of is that it depends on 
the spin polarization Ad. For 7^ K^{uo) is finite 
at all values of uo and the collision integral is well de- 
fined. A peculiar situation arises when A// = 0. In 



this case K^{u) cx 



is singular at a; = 0. At 



the first glance there is nothing special in such a de- 
generacy which is the same as for the standard case 
of the electron-electron interaction in the singlet chan- 
nel^'^. The difference, however, is in the form of I^^ 



as compared to Tc 



(1,2) 



As is seen from Eqs. (|4])-(|6]), 



the cancellations of the "in" and "out" terms lead to: 



Fa ' = 0) = for any /^(F), while = 0) 7^ 

unless /dif = /t(F) — /|(F) is identically zero. This 
means that the singularity Kp{(jo) oc ^jp72 is not danger- 

ous for terms proportional io la (which are the only 
terms that arise in the case of electron-electron inter- 
action in the singlet channel) but it leads to the diver- 

f3) 

gency of the term cx la corresponding to the spin-fiip 
processes due to the electron-electron interaction in the 
triplet channel. So we have an infrared catastrophe in the 
collision integral in the case of a spin-dependent electron 
energy distribution with A/1 = 0. 

5. Derivation of the collision integral. Before we 
discuss the physical origin of such a catastrophe and the 
ways to cure it we briefly outline the derivation of the 
collision integral Eqs. (|4|)-(|6l).P- (fTQ|) . 

Following the original work of Keldysbi^ we represent 
the collision integral as follows: 



^coll = —i{^G — G^)l2 



where 







g^ 



(18) 



(19) 
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and S^'^'^ and Q^^^^^ are retarded, advanced and 
Keldysh components of the self-energy part and an exact 
single-particle Green's function. In the two-quasiparticle 
collision approximation we adopt in this paper, the self- 
energy part is given by: 





= I>J 






(20) 




= 


^E+u - 




(21) 




= 


^E+oj ~ 




-Gi+J(22) 



where the integration over co and a proper summation 
over spin indices a, 7, 6 is assumed in all three equa- 
tions. 

In this approximation we neglect (a) the interaction 
corrections to the vertex part F = 1 and (b) the inter- 
action corrections to the single-particle Green's function 
Ge-\-uj which is supposed to be equal to the corresponding 
Green's function without interaction Ge+uj- In Eqs. ([2Q|) 
we denote by the dynamically screened interaction: 



(23) 



where a* {i = 1, 2, 3) are the Pauli matrices d-x^y^z and 
is the unit matrix. 

The 4x4 interaction matrix Df^^^^ obeys an RPA-like 
equation: 



D, 



R{A) 



UAj^UAkI^ll^^Df^^\ (24) 



where Ui {i = 1,2,3) is the bare interaction constant 
= F/ {4:u) of the electron-electron interaction in the 
triplet channel: 



(25) 



Pi,2jq<q* j=l,2,3 af3^S 

and U4 = = F^ -\- V{q) is the bare interaction in 
the singlet channel which at small q is dominated by the 
Coulomb interaction V{q) 00: 



Pi,2,q<q* ^7 



^[V{q)^F^] (26) 
X ii^iiPi) ^a{pi - q)] [i^\{P2) ^j{p2 + q)]' 



The generalized polarization bubbles are given by 



n 



R,A,K 
kl 



a,f3 



(27) 



where the retarded and advanced polarization bubbles 
~^ ^aa^^^J/3+^aa^i^/3^^ for the spin-independent 
single particle Hamiltonian (no magnetic impurities and 
no spin-orbit interaction) contain retarded, advanced and 
Keldysh Green's functions. 



The Keldysh component Df^ of the dynamically 



R 



screened interaction is expressed explicitly through Dij 



and Dfji 



(28) 



where 11^ is given by Eq. ([27|) but with tt^ 



r<R nA 



r<A r'R 



The central point of the derivation of the collision inte- 
gral is the ansatz that involves the non-equilibrium elec- 
tron energy distribution function hcr^E- 



G 



K 

E,a 



(G 



R 

E,a 



^E,a) ^cr,E' 



(29) 



The similar ansatz has been suggested by Keldysh^^ for 
the exact Green's functions. We will be using Eq. ([29|) 
for the Green's functions without electron- electron inter- 
action. The reason is that the perturbation theory in 
interaction can be built using Eq. ([29|) with an arbitrary 

"initial" distribution function h^^"^^ compatible with the 
Fermi statistics. It will cancel out anyway in the final 
result as the initial distribution has to be forgotten in 
the non-equilibrium steady state. Diagrammatically this 
cancellation happens because of the proliferation of sin- 
gular "loose diffusons"^^ which, however, make impossi- 
ble the perturbative analysis. There is only one single 
choice of h^^^^ in Eq. ([29|) - the true solution h(j,E of the 
kinetic equation - when such proliferation does not oc- 
cur and all the diagrams with loose diffusons are equal 
to zeroi^. 

Using Eq. ([29|) one can obtain the following expressions 
for the polarization bubbles: 



-i [{ha^E'+u. - hp.E')G%^^Gi. - 2iv\ ,(30) 
M {(hc^E'^u. - hp.E')G%+^G% + 2iv] ,(31) 

i (1 - hoc,E'+u;h(3^E')^GE'+u;^GE', (32) 



where /S.Ge = {G^ — G^)e and an integration over 
J dF' is assumed. The next step is the standard dis- 
order average of the product G^,j^^G^, with the result 



i'+a.Gi)c.,q = 27n./(Z)q2-iu;). 



(33) 



This eliminates the dependence on F' everywhere but in 
the distribution functions ha,E' • 

At this point it is appropriate to note on the difference 
between iTaa^^ (or the spin-independent case) and tt^^"^^ 
with a /3. In the former case one can use an identity 

J dF' {K^E'+o. - K^E') = 2cJ, (34) 

which holds for an arbitrary function ha,E' which at 
\F'\ oo converges sufficiently fast to sign(£^'). Then 
one immediately obtains the standard result which is in- 
dependent of the electron energy distribution function: 



<. = (^^)*. (35) 
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In contrast to that, the corresponding integral in 7r^_Q, 
contains a part / dE' {ha,E' —h-a,E') proportional to the 
total spin polarization Ai of the non-equilibrium state 
given by Eq. ([3]). Thus we obtain: 



vM 



L>q2 



(36) 



where the sign ^ corresponds to a =] {[). For complete- 

K . 



ness we also give an expression for tt^: 



-2iv 



(1-/1, 



(37) 

The dependence of Eqs. ([36|) . ([37|) on the spin projections 
Of, p makes the 4x4 matrices n^^'"^'^ non-diagonal: 



n 



R,A,K 



( n2 

^3 



V 



-zEs 

n2 

no Hi 

Hi no 



R,A,K 



where for any of the omitted superscripts A, K 



Ho 

n2 



Hi 



(38) 



(39) 



Correspondingly, the 4x4 matrices Df-^'^ found from 
Eqs. ([24 |) . ([28|) appear to have the off-diagonal structure 
similar to Eq. For D^^^^ = (1 - V^^^^^Y^V we 

obtain: 



D 



R,A,K 





D2 


-iD3 








iD3 


D2 


















Dzo 


\ 













R,A,K 



(40) 



K 



D 



^00 



(l-2F7r^)(l-2F7r^)' 
{l-2UP7T^){l-2UP7r^) ' 47r^7r^ 



hi 



D 



K ^ 

(l-2F7r^)(l-2/7^7r^) 



Caf3 



(l-2/7^7r^)(l-2F7r^) 
1 



27r^(l-2F7r^)' 

27r^(l-2F7r^)' 
(43) 



In Eqs. (|4lj),(|43|) we took the limit ^ oo that cor- 
responds to |V(q)7r^*^^^(cc;, q)| 1 which is always the 
case at small enough |q|, as V(q 0) diverges because 
of the long-range character of Coulomb interaction. 

Substituting Eqs. ([20])-(|43l) into Eq. ([18]), neglecting 
interaction correction to the Green's function Qe ^ Ge 
and using Eq. (|33|) for the disorder average of still not 
averaged pair of Green's functions we arrive at the main 
quantitative result of this paper given by Eq. (jlj)-([TO|). 

6. Conclusion. We have shown above that the relax- 
ation of the spin-dependent electron energy distribution 
at the total spin magnetization = differs qualita- 
tively from the case 7^ 0. It is only in this case that 
the complete relaxation to a spin-independent Fermi dis- 
tribution is possible due to electron-electron interaction 
alone. And it is in this case that the infrared catastro- 
phe is encountered in the collision integral for a quasi- Id 
disordered wire. As a result, an anomalously fast relax- 
ation to a spin-independent non-equilibrium distribution 
happens well before the complete equilibrium is reached. 
The corresponding collision integral responsible for such 
a fast relaxation can be approximated as 



with 



R{A) 
2 

R{A) 



F^(l 



R{A] 



(41) 



^R{A) 



1 



(l-2F-^f/^^)(l 



2F^7r 



R{A). 

it ^ 



D 



RiA) 
00 — 

R(A) 
zO — 



1 - 2F^7r^(^) ' 
UP 

1 - 2/7^7r^(^) " 



(42) 



1 



27r^(^) ' 



RiA) 



0. 



Respectively, for D we have: 

^3 = (^nCn -^^TCiT)^^ 



(44) 



where 




jil- K,E'h-^,E')dE' , (45) 



with C{F) = 

A remarkable feature of Eq. (|45|) emerging due to the 
infrared catastrophe is the quasi-elastic form of the col- 
lision integral. If the infrared cut-off 1/ro in the diver- 
gent integral Ji^^^ \oo\^/'^ small compared to the effec- 
tive temperature, one can set cc; = in the distribution 
functions fa- entering Eq. (|6]). Thus we arrive at Eq. (|44|) 
which dependence on fa{E) is identical to the elastic part 
of the collision integral due to magnetic impurities. The 
only difference is that the coefficient i^cr also depends on 
the integral of the distribution functions. Thus the triplet 
part of the electron-electron interaction acts in this case 
similar to the magnetic impurities. 
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